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Statistical Significance of the Gallium Anomaly 
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We calculate the statistical significance of the anomalous deficit of electron neutrinos measured in 
the radioactive source experiments of the GALLEX and SAGE solar neutrino detectors taking into 
account the uncertainty of the detection cross section. We found that the statistical significance of 
the anomaly is about 3.0a. A fit of the data in terms of neutrino oscillations favors at about 2. To- 
short-baseline electron neutrino disappearance with respect to the null hypothesis of no oscillations. 



PACS numbers: 14.60.Pq, 14.60.Lm, 14.60.St 

The GALLEX [J-il and SAGE d-Q Gallium solar 
neutrino experiments have been tested with intense ar- 
tificial ^^Cr and '^'^Ar radioactive sources placed inside 
the detectors. The results of these "Gallium radioactive 
source experiments" indicate a ratio R of measured and 
predicted ''^Ge event rates which is smaller than unity: 



= 0.953 ±0.11 
i?g2 = 0.812l°;J?, 
= 0.95 ±0.12, 

pS2 _ n 7QT +0.084 
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where Gl and G2 denote the two GALLEX experiments 
with ^^Cr sources, SI denotes the SAGE experiment with 
a '"^^Cr source, and S2 denotes the SAGE experiment with 



37 



Ar source. 



Assuming Gaussian probability distributions and tak- 
ing into account the asymmetric uncertainties of and 
we have the probability distributions shown by the 
dashed, dotted, dash-dotted and dash-dot-dotted lines in 
Fig. [T] The combined probability distribution PuGb-It) 
shown in Fig. [T] gives the average ratio 



pGa _ n nf^+0. 05+0. 10+0. 15 
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(5) 



where the uncertainties are at 68.27% C.L. (Ict), 95.45% 
C.L. {2a), 99.73% C.L. (3cr). Thus, the number of mea- 
sured events is about 2.8a smaller than the prediction. 
This is the "Gallium anomaly" . 

As indicated by the "B" subscript, the ratios in 
Eqs. (HI)-® have been calculated with respect to the rate 
estimated using the best-fit values of the cross section of 
the detection process 
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calculated by Bahcall 



a^TCr) = {58.1111) 
a^^(37Ar)= {TO.Otti) xlO-^Scm^ 



X 10"''^ cm^ , 



(7) 
(8) 



The uncertainties of these cross sections are not taken 
into account in the experimental ratios in Eqs. (II|)-(|4l). 
These uncertainties are large because only the 

cross section of the transition from the ground state of 
^^Ga to the ground state of ''"'^Ge is known with preci- 
sion from the measured rate of electron capture decay of 
^^Ge to ^^Ga. Electron neutrinos produced by ^^Cr and 
■^''Ar radioactive sources can be absorbed also through 
transitions from the ground state of ^^Ga to two ex- 
cited states of ^^Ge at 175 keV and 500 keV, with cross 
sections which are inferred using a nuclear model from 
p + ""^Ga ^^Ge + n measurements [Tlj |. 

Hence, at least part of the deficit of measured events 
with respect to the prediction could be explained by an 
overestimation of the transitions to the two excited states 
office ia [T^. However, since the contribution of the 
transitions to the two excited states of ''^Ge is only 5% 
[1], even the complete absence of such transitions would 
reduce the ratio of measured and predicted ^^Ge event 
rates to about 0.91 ± 0.05, leaving an anomaly of about 
1.7a [13. 

We think that for a correct assessment of the statistical 
significance of the Gallium anomaly simple approaches 
based on either accepting the Bahcall cross section in 
Eq. ([7]) without taking into account its uncertainty or 
suppressing without theoretical motivations the transi- 
tions to the two excited states of ^^Ge are insufficient. 
A correct assessment of the statistical significance of the 
Gallium anomaly can be done by taking into account 
the large uncertainties of the transitions to the two ex- 
cited states of ""^Ge (SrrlOj. The most reliable estimate of 
these transitions and their uncertainties have been done 
by Haxton in Ref. [l3|, leading to the total cross section 
for a ^^Cr source 



crH(^'Cr) = (63.9 ±6.8) x 10 



-46cm2. 



(9) 
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Notice that the average value of this cross section is even 
larger than the Bahcall cross section in Eq. ([7]). This 
leads to an enhancement of the Gallium anomaly. How- 
ever, the uncertainty of ani^^Cr) is rather large. Hence, 
a correct assessment of the statistical significance of the 
Gallium anomaly requires an accurate treatment of the 
cross section uncertainty. 

Since the ratios in Eqs. (H])-© have been calculated 
with respect to the best-fit value in Eq. ([7]) of the Bah- 
call cross section for a ^^Cr source, these ratios must be 
rescaled by 



^H(^^Cr) 
<(5iCr) 



1.10±0.12. 



(10) 



For the SAGE ■^''Ar source experiment, we evaluate 
the detection cross section and its uncertainty as follows. 
The cross section is given by Q 



cr(^^ Ar) = cTgJ l-f 0.695 



BGTi 



75 



BGT 



0.263 



BGT 



500 



BGT. 



(11) 

where dgs = 66.2 x 10 '^^ cm^ is the cross section from the 
ground state of ^^Ga to the ground state of ""^Gc, BGTgg 
is the corresponding Gamow- Teller strength and BGT175 
and BGT500 are the Gamow- Teller strengths of the tran- 
sitions from the ground state of ''"'^Ga to the two excited 
states of ^^Ge at 175 keV and 500 keV. The coefiicients 
of BGTi75/BGTgs and BGTgoo/BGTg^ are determined 



by phase space. In Ref. [l^, Haxton estimatecfl 

BGTi75/BGTg, = 0.19 ± 0.18 , 
BGTsoo/BGTgs = 0.13 ± 0.02 . 



Thus, we obtain 



(77.3 ±8.2) X 10 



-46cm2, 



and 



i?g(3^Ar) 



^H (^'Ar) 
(37Ar) 



1.10 ±0.12. 



(12) 
(13) 
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which has the same value of i?g(5iCr) in Eq. There- 
fore, all the ratios in Eqs. ©-(Ill) must be rescaled by 
i?H = i?H(51Cj.) = i?H(37Ar). 

One must also take into account that the value of the 
cross section is bounded from below by the cross section 
CTgs of the transition from the ground state of ^^Ga to the 
ground state of ^^Ge 



^bf 



0.95. 



(16) 




FIG. 1. Solid lines: probability distributions pjjGa(r) and 

P/{Ga(r), as indicated by the labels. Dashed, dotted, dash- 
dotted and dash-dot-dotted lines: probability distributions 
PHgi(^). Pflg2(?-), p^jsi(r) and P^jS2(r), respectively. 



In the following we calculate the probability distribu- 
tion of 



(17) 



by taking into account the uncertainty of the denomina- 
tor i?g given in Eqs. (|10|) and (|15p . This is the theoretical 
uncertainty of the cross section that has not been taken 
into account in the ratios ([T])-©, which have been evalu- 
ated using the best- fit values of the Bahcall cross sections 
in Eqs. (O and ([8|). We assume a Gaussian probability 
distribution truncated below i?g : 



exp 
0, 



Ai?H 



r>R 



B ' 



(18) 



r < R 



B ' 



with = 1.10 and Ai?g = 0.12. 

The probability distribution of the ratio i?*^** 
Eq. ([HI) is given by (see Section 2.4.4 of Ref. fli |) 



PflGa(r) = / p^Ga(rs)pflH(s)sds. (19) 
JrS 



In Ref. the values of BGTirs/BGTgs and BGTsoo/BGTgs 
can be extracted, respectively, from Eqs. (12) and (7), taking 
into account Eq. (1). As explained by Haxton, BGT175 requires 
a theoretical calculation, whereas for BGT500 it is reasonable to 
adopt the corresponding (p, n) value. 



Figure [T] shows the probability distribution PnOn (r) of 
R^^ derived from the experimental data in Eqs. (IT])-(|3]) 
and the result of the integral in Eq. p^ . One can see 
that p7jGa(r) is peaked at a smaller value than p^Ga(r'), 
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but the uncertainty is larger. We obtain 

nGa _ n 7f5+0. 09+0. 17+0.24 
n — U. (0_Q 08-0.15-0.21 



(20) 



where the uncertainties are at 68.27% C.L. (la), 95.45% 
C.L. {2a), 99.73% C.L. (3ct). From a comparison of 
these uncertainties and from Fig. [T] one can see that 
the probabihty distribution is approximately Gaussian, 
with slightly asymmetric uncertainties and tails which 
decrease slightly faster than Gaussian tails. 

The probability of R^" < 1 is 99.86% (3.0ct anomaly), 
slightly larger than the probability of R^^ < 1 , which is 
99.75% (2.8cr anomaly). Therefore, the Gallium anomaly 
remains statistically significant after taking properly into 
account the cross section uncertainty. 

For the four individual Gallium radioactive source ex- 
periments, using the same method as above, from the 
experimental values in Eqs. (HJ-Q we obtain 





nGl / pH 


n Q/l+0. 13+0. 26+0. 40 
— U.04:_Q j^2-0.23-0.33 ' 


(21) 




nG2 / nH 


r, 7-1 +0.12+0.24+0.36 
~ ^- ' ^-0.11-0.21-0.31 ' 


(22) 




pSl /pH 


p. o/i +0.14+0.28+0.42 
— U.O'±_o 23-0.24-0.35 ' 


(23) 




pS2 / pH 

— /-Kb 


r, yrv + 0. 10+0. 21+0. 31 

— <U_o 09-0.17-0.25 ' 


(24) 



with IfT, 2iT, 3(7 uncertainties. A comparison of these un- 
certainties shows that the probability distributions are 
approximately Gaussian, with slightly asymmetric un- 
certainties. 

Since the Gallium anomaly is confirmed by the new 
statistical analysis which takes into account the uncer- 
tainty of the detection cross section, it is plausible that it 
is due to a physical mechanism. In the following, we con- 
sider the possibility of electron neutrino disappearance 
due to short-baseline oscillations [l^ [l5l - l2l| (another ex- 
planation based on quantum decoherence in neutrino os- 
cillations has been proposed in Ref. p^). 

We consider the electron neutrino survival probability 



P^^^^^ {L, E) = l- sin^ 2i9 sin" 



Am'^L 

4:E 



(25) 



where 'd is the mixing angle, Am^ is the squared- 
mass difference, L is the neutrino path length and E 
is the neutrino energy. This survival probability is effec- 
tive in short-baseline (SBL) experiments in the frame- 
work of four-neutrino mixing schemes (see Refs. psl - 
[2^), which are the simplest extensions of three- neutrino 
mixing schemes which can accommodate the two mea- 
sured small solar and atmospheric squared-mass differ- 
ences Am§QL ^ 8 X lO"'^ eV^ and Ato^^m - 2 x 10"^ eV^ 
and one larger squared-mass difference for short-baseline 
neutrino oscillations, Am^ ^0.1 eV^. The existence of a 
fourth massive neutrino corresponds, in the flavor basis, 
to the existence of a sterile neutrino Vg. 

We performed a maximum likelihood analysis (see 
Ref. [231 ) of the Gallium data as foUowfl We started 



with the calculation, for each experiment, of the value of 
the ratio of the event rate as a function of sin^ 2i? and 
Am^ and the event rate in absence of neutrino oscilla- 
tions (see Ref. for details): 



R'^isin'^ 2i^, Am^) = 



kdVL~^E.b'[a'[P^^J;,^{L,E.) 



(26) 

where the index k labels the experiments (fc = 
Gl, G2, SI, S2), the index i labels the lines emitted 
in ^^Cr or ^^Ar electron captures with energies Ei, 
and erf are the correspondin g b ranching ratios and cross 
sections (see Table I of Ref. [13): L is the neutrino path 
length and J^, dV is the integral over the volume of each 
detector (see Table II of Ref. (H). 

The uncertainty of i?^ is correlated in the calculation 
of the combined probability distribution of the four ex- 
perimental ratios in Eqs. (I?n) - ([M|) . Using a method sim- 
ilar to that utilized for the derivation of Eq. ([TO)) (see 
Section 2.4.4 of Ref. [l3l), we obtain the combined prob- 
ability distribution 



B L fe 



Pi,^{s)s^ds, (27) 



where R = {R'^\ R'^^ , R^\ R^^) and f = 
(r^\r^2 ,.31^^32^^ The authors of Ref. [28] con- 
sidered a correlation of the systematic errors of the two 
GALLEX experiments and the two SAGE experiments. 
Since such correlation is not documented in the ex- 
perimental publications, where the combined ratio was 
calculated as a weighted average, without correlations, 
we adopt the conservative approach of considering the 
systematic experimental errors as independenllf|. 

The likelihood function of the oscillation parameters 
sin^ 2i^ and Am^ is given by 



Cisi-a^ 2d, Ato^) = p^(i?(sin2 2z9, Am^)) , 



(28) 



with the four components of i?(sin^ 2i9, Am^) given by 
Eq. (|26l) . Figure [H shows the allowed regions in the 
sin^ 2??-Am^ plane and the marginal Ax^ = ~ Xmin'^ 
for sin^ 2i9 and Am^, from which one can infer the cor- 
responding uncorrelated allowed intervals. In the max- 
imum likelihood analysis x^(sin^ 2i9, Am^) is given by 
-21n/:(sin2 2'd, Am^) + constant. 

The best-fit values of the oscillation parameters are 

(29) 



sin" 2^9 



bf 



0.50, 



Ami,- 



2.24 eV 



A standard least-squares analysis would lead to misleading re- 



sults, because it does not allow us to take into account the lower 
bound in Eq. ^ for 
^ A correlation of the systematic experimental errors can be taken 
into account in Eq. II27I I by replacing YlkPR.'' {^*'') with a mul- 
tivariate Gaussian distribution with the appropriate covariance 
matrix. 
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Gallium 

68.27% C.L (1a) 
90.00% C.L. 
95.45% C.L. (2a) 

1.00% C.L. 

1.73% C.L. (3a) 




of sin^ 2d and Am^ indicate thal|3 



sin 2-d 



FIG. 2. Allowed regions in the sin^ 2i9-Am^ plane and 
marginal Ax^'s for sin^ 2i? and Am^ obtained from the com- 
bined fit of the results of the two GALLEX ^^Cr radioactive 
source experiments and the SAGE ^^Cr and '^^Ar radioactive 
source experiments. The best-fit point corresponding to Xmin 
is indicated by a cross. 



The value of the likelihood ratio between the null hy- 
pothesis of no oscillations and the oscillation hypothesis, 



£0 



£(sin^ 2i9bf, Am^^) 



= 8 X 10" 



(30) 



is in favor of the oscillation hypothesis. It corresponds 
to Ax^ = 9.7, which, with two degrees of freedom, dis- 
favors the null hypothesis of no oscillations at 99.23% 
C.L. {2.7a). The small difference between this statisti- 
cal significance of the indication in favor of the Gallium 
anomaly and that obtained from Eq. (1^0]) (3. Oct) is due to 
the different analysis of the data. Although the neutrino 
oscillation analysis leads to a better fit of the four data 
in Eqs. (PT|) -([M 1) (the best-fit values of the oscillation 
parameters in Eq. ^ give = 0.75, i?^^ = 0.75, 
i?si ^ 0.73 and R^"^ = 0.72), the correlation of the theo- 
retical uncertainty of iJg slightly disfavors a fit in which 
the deviations of the data from the best-fit values do not 
have the same sign. 

From Fig.[5]one can see that the marginal distributions 



sin^ 2d > 0.07 , 



Ato^ > 0.35 eV^ 



(31) 



at 99% C.L.. These bounds indicate that the short- 
baseline disappearance of electron neutrinos may be 
larger than that of electron antincutrinos, which is 
bounded by the results of reactor neutrino experiments 
p, 19, 20]. This could be an indication of a violation 
of the CPT symmetry ^ (CPT implies that P^^^u^ = 
Pu^-^vc, for any flavor a = e,ii, r; see Ref. [l^l). However, 
according to a recent calculation [3l| the De fluxes pro- 
duced in nuclear reactors are about 3% larger than the 
standard ones used in the analysis of reactor antineutrino 
data (see Ref. [l^l). A comparison of the new reactor an- 
tineutrino fluxes with the data of several reactor neutrino 
experiments suggests the existence of a reactor antineu- 
trino anomaly [2gl which is compatible with the Gallium 
anomaly in a standard CPT-invariant neutrino oscilla- 
tion framework. In this case, the indication in favor of 
CPT violation obtained by comparing the results of the 
neutrino oscillation analysis of Gallium and reactor data 
is weakened, but the plausibility of the existence of a Gal- 
lium anomaly is reinforced by its compatibility with the 
reactor antineutrino anomaly. 

CPT violation in short-baseline electron neutrino dis- 
appearance can be tested with high accuracy in future ex- 
periments with pure and well-known i^e and Ve beams, as 
beta-beam [33] and neutrino factory [U,!!!] experiments. 
Although the possibility of CPT violation is theoretically 
problematic [36[, it cannot be dismissed in phenomeno- 
logical analyses of experimental results. It is interesting 
to notice that recently another indication of a violation of 
the CPT symmetry has been found in the MINOS long- 
baseline i^fj and Df^ disappearance experiment (stI . Issj . 

There is also a growing experimental interest in favor of 
possible tests of the Gallium anomaly. In addition to the 
future experimental possibilities to test the short-baseline 
disappearance of electron neutrinos discussed in Ref. [l3j , 
the authors of Ref. [2l| presented recently a plan to make 
an improved direct measurement of the Gallium anomaly 
with the liquid Gallium metal used in the SAGE exper- 
iment and a new vessel divided in two zones, which can 
measure a variation of the electron neutrino disappear- 
ance with distance. The Borexino collaboration is study- 
ing the possibility of a radioactive source experiment j39[ 
which could provide a "smoking gun" signal by measuring 
the oscillation pattern inside the detector. Other possi- 
ble measurements with radioactive sources and different 
detector types has been recently discussed in Ref. 40]. 

The existence of at least four massive neutrinos, one of 
which has a mass larger than about 0.6 eV in order to gen- 
erate the squared-mass difference in Eq. (pij) . can have 



These bounds are weaker than those presented in a previous 
version of this paper (arXiv: 1006.3244v2) in which the corre- 
lation of the uncertainty of i?g in the calculation of the com- 
bined probability distribution of the four experimental ratios in 
Eqs. II21I I- II24I I was not taken into account. 
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important implications for cosmology (see Refs. |4ll443| ). 
The current indications of cosmological data analyzed in 
the framework of the standard cosmological model are 
controversial. On one hand, there are indications that 
the effective number of neutrino species r nay be larger 
than three from Big Bang Nucleosynthesis [43| and from 
the Cosmic Microwave Background Radiation [4^. This 
is consistent with a thermalization of sterile neutrinos 
due to active-sterile oscillations before Big Bang Nu- 
cleosynthesis induced by the large values of the mixing 
parameters in Eq. ([5T]) (46| . On the other hand, anal- 
yses of Cosmic Microwave Background Radiation data 
and Large Scale Structure data constrain the mass of 
a fourth thermalized neutrino to be smaller than about 
0.7 eV j47l - l50l |. Hence, either the heavy neutrino mass is 
close to the standard cosmological bound or the existence 
of short-baseline neutrino oscillations is connected with 
non-standard cosmological effects, as those discussed in 



Refs. Hil-il. 

In conclusion, we have estimated the uncertainty of the 
deficit of electron neutrinos measured in the radioactive 
source experiments of the GALLEX and SAGE 0- 
0| solar neutrino detectors taking into account the uncer- 
tainty of the detection cross section estimated by Haxton 
in Ref. fld\. The result shows that the Gallium anomaly 
is statistically significant, at a level of about 3.0a. The 
analysis of the data in terms of neutrino oscillations in- 
dicates values of the oscillation amplitude sin^ 2d > 0.07 
and squared-mass difference Ato^ > 0.35 eV^ at 99% 
C.L.. 
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